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We present a gauge-invariant formalism to study the evolution of the curvature and entropy 
perturbations in the case in which spatial and time variations of the inflaton decay rate 
into ordinary matter are present. During the reheating stage after inflation curvature 
perturbations can vary with time on super-horizon scales sourced by a a gauge-invariant 
inflaton decay rate perturbation. We show that the latter is a function not only of the spatial 
variations of the decay rate generated during inflation, as envisaged in a recently proposed 
scenario, but also of the time variation of the inflaton decay rate during reheating. If only 
the second source is present, the final curvature perturbation at the end of the reheating 
stage is proportional to the curvature perturbation at the beginning of reheating with a 
coefficient of proportionality which can be either smaller or larger than unity depending 
upon the underlying physics governing the time variation of the inflaton decay rate. As 
a consequence, we show that the standard consistency relation between the amplitude of 
curvature perturbations, the amplitude of tensor perturbations and the tensor spectral 
index of one-single field models of inflation is violated and there is the possibility that the 
tensor-to-curvature amplitude ratio is larger than in the standard case. 
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1 Introduction 



One of the basic ideas of modern cosmology is that there was an epoch early in the history 
of the universe when potential, or vacuum, energy associated to a scalar field, the inflaton 
0, dominated other forms of energy density such as matter or radiation. During such 
a vacuum-dominated era the scale factor grew exponentially (or nearly exponentially) in 
time. During this phase, dubbed inflation pflE], a sixxa.ll , smooth spatial region of size less 
than the Hubble radius could grow so large as to easily encompass the comoving volume of 
the entire presently observable universe. If the universe underwent such a period of rapid 
expansion, one can understand why the observed universe is so homogeneous and isotropic 
to high accuracy. 

Inflation has also become the dominant paradigm for understanding the initial condi- 
tions for structure formation and for Cosmic Microwave Background (CMB) anisotropy. In 
the inflationary picture, primordial density and gravity-wave fluctuations are created from 
quantum fluctuations "redshifted" out of the horizon during an early period of superlumi- 
nal expansion of the universe, where they are "frozen" 012130 El- Perturbations at the 
surface of last scattering are observable as temperature anisotropy in the CMB which was 
first detected by the Cosmic Background Explorer (COBE) satellite [HI El EH] • The last and 
most impressive confirmation of the inflationary paradigm has been recently provided by 
the data of the Wilkinson Microwave Anistropy Probe (WMAP) mission which has marked 
the beginning of the precision era of the CMB measurements in space jllj . The WMAP 
collaboration has produced a full-sky map of the angular variations of the CMB, with un- 
precedented accuracy. WMAP data confirm the inflationary mechanism as responsible for 
the generation of curvature (adiabatic) super-horizon fluctuations. 

Despite the simplicity of the inflationary paradigm, the mechanism by which cosmo- 
logical adiabatic perturbations are generated is not yet fully established. In the standard 
picture, the observed density perturbations are due to fluctuations of the inflaton field it- 
self. When inflation ends, the inflaton oscillates about the minimum of its potential and 
decays, thereby reheating the universe. As a result of the fluctuations each region of the 
universe goes through the same history but at slightly different times. The final temper- 
ature anisotropies are caused by the fact that inflation lasts different amounts of time in 
different regions of the universe leading to adiabatic perturbations. Under this hypothesis, 
the WMAP dataset already allows to extract the parameters relevant for distinguishing 
among single- field inflation models [T2]. 
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An alternative to the standard scenario is represented by the curvaton mechanism [T3[ 
ITU ITo] where the final curvature perturbations are produced from an initial isocurvature 
perturbation associated to the quantum fluctuations of a light scalar field (other than the 
inflaton), the curvaton, whose energy density is negligible during inflation. The curvaton 
isocurvature perturbations are transformed into adiabatic ones when the curvaton decays 
into radiation much after the end of inflation. 

Recently, another mechanism for the generation of cosmological perturbations has been 
proposed [TBI [XH] ■ It acts during the reheating stage after inflation and it was dubbed the 
"inhomogeneous reheating" mechanism in Ref . [IB] . The coupling of the inflaton to normal 
matter may be determined by the vacuum expectation value of fields x's of the underlying 
theory. If those fields are light during inflation, fluctuations Sx ~ where H is th Hubble 
rate during inflation, are left imprinted on super-horizon scales. These perturbations lead 
to spatial fluctuations in the decay rate T of the inflaton field to ordinary matter, ~ 

causing adiabatic perturbations in the final reheating temperature in different regions of 
the universe. 

Interestingly, these different scenarios have different observational predictions. The 
curvaton scenario and the one based on variation of the decay rate allows to generate 
the observed level of density perturbations with a much lower scale of inflation and thus 
generically predicts a smaller level of gravitational waves. Furthermore, because the field 
responsible for the fluctuations is not the inflaton, it can have significantly larger self 
couplings and thus density perturbations could be non-Gaussian. The non-Gaussianity 
can be large enough to be detectable by CMB and Large Scale Structure observations 
contrary to what predicted in the traditional one-single field model of inflation, where the 
level of non-Gaussianity is very small [T§] . 

Contrary to the standard picture, both the curvaton and the inhomogeneous reheat- 
ing mechanism exploit the fact that the total curvature perturbation (on uniform density 
hypersurfaces) ( can change on arbitrarily large scales due to a non-adiabatic pressure 
perturbation which may be present in a multi-fluid system |2Til I2*T| [221 122] ■ While the en- 
tropy perturbations evolve independently of the curvature perturbation on large scales, the 
evolution of the large-scale curvature is sourced by entropy perturbations. 

In this paper, we present a gauge- invariant formalism to study the evolution of the gauge- 
invariant curvature and entropy perturbations in the inhomogeneous reheating mechanism. 
To do so, we extend the gauge- invariant formalism first introduced in Ref. [21j, which is 
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appropriate to describe the evolution of curvature and entropy perturbations in multi-fluid 
cosmologies when energy transfer between fluids is included. We show that the curvature 
perturbation during the reheating stage can vary with time on super-horizon scales sourced 
by a a gauge-invariant inflaton decay rate perturbation 

ST G1 = ST-t 5 -^, (1.1) 

where p^ is the inflaton energy density and Sp^ its perturbation. The gauge-invariant 
inflaton decay rate perturbation is a function not only of the spatial variation of the decay 
rate ST generated during inflation, but also of the time variation of the inflaton decay rate 

T. To our knowledge, this new source proportional to T has never been discussed in the 
literature. 

As an application, we study the evolution of the curvature perturbation in the inhomo- 
geneous reheating scenario where fluctuations of the inflaton decay rate ST are induced by 
the fluctuations of some light scalar field during inflation and confirm the results of Refs. 

We extend their findings showing that variation of the total curvature perturbation 
( on super-horizon scales take place even when the inflaton decay rate changes with time 
during reheating. This new effect is proportional to the total curvature perturbation gen- 
erated during the inflationary period ( in and causes either an increase or a depletion of ( in 
depending upon the underlying model responsible for the time variation of T. Furthermore, 
we show that the standard consistency relation between the amplitude of curvature per- 
turbations, the amplitude of tensor perturbations and the tensor spectral index is violated. 
Finally, for completeness, we extend our analysis to the case in which the classical inflaton 
field does not release its energy perturbatively, but very rapidly (explosively) decays into 
either its own quanta or into other bosons due to broad parametric resonance, the so-called 
preheating stage. 

2 The perturbative inflaton decay 

Our starting point are the equations governing the evolution of cosmological perturbations 
during the reheating stage. We follow the gauge-invariant approach developed in Ref. 
for the general case of an arbitrary number of interacting fluids in general relativity. 

At the end of inflation, once the Hubble rate drops below the mass of the inflaton field 
</>, the inflaton starts oscillating around the minimum of its potential. Averaged over several 
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oscillations, the effective equation of state is (P^/p^) = 0, where P^ and p^ are the inflaton 
pressure and energy density, respectively. The coherent oscillations of the inflaton field are 
equivalent to a fluid of non-relativistic particles [22] • The vacuum energy during inflation 
is transformed into the energy density of the coherent inflaton oscillations. Assuming the 
inflaton is unstable and decays into light particles ("radiation") with a decay rate T, this 
represents an energy transfer from the pressureless inflaton fluid to the radiation fluid with 
energy density p 7 . 

The evolution of the background FRW universe during the reheating stage is governed 
by the Friedmann constraint 

& = "fp, (2-1) 
H = -AttG (p + P) , (2.2) 

and the continuity equation 

p = -3H(p + P), (2.3) 

where the dot denotes differentiation with respect to the coordinate time t, H = a/ a is the 
Hubble parameter, and p and P are the total energy density and the total pressure of the 
system. The total energy density and the total pressure are related to the energy density 
and pressure of the inflaton field and radiation by 



P = P^ + Pj, 

P = P + P 7 , (2.4) 

where P 7 is the radiation pressure. The inflaton field and the radiation component have 
energy- momentum tensor Tj*K and T^, respectively. The total energy momentum tensor 
rpfiu _ rppv _|_ rpLw j g covariantly conserved, but we allow for energy transfer between the 
fluids, 

= Qfa, (2.5) 

where Qfy and Q^a are the generic energy-momentum transfer to the inflaton and radiation 
sector respectively and are subject to the constraint 

%) + %) = • (2-6) 
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The continuity equations for the energy density of the inflaton field p^ and radiation p 7 in 
the background is thus [2H] {Q<j> = Q ^, Q 7 = Q° 7) ) 



p^ = -3H (p + P$) + , 

p 7 = -3H (p 7 + P 7 ) + Q 7 . (2.7) 

From now on, we parametrize, the energy transfer from the inflaton to radiation by [21] 

Q<t> — —Fp<f> , 

Q 7 = r P0 , (2.8) 

where T is the decay rate of the inflaton into radiation. The positions (|2.8|) are valid in the 
case in which the inflaton decays into light states through a perturbative process. However, 
if the classical inflaton field very rapidly (explosively) decays into either its own quanta 
or into other bosons due to broad parametric resonance, the so-called preheating stage [2H] , 
Eqs. ()2.8|) should be modified. The corresponding equations will be discussed in Section |3] 

The energy conservation equations are therefore 

P4> = - P<t} (3H + T) , (2.9) 

p 7 = -4#p 7 + r/v, (2.10) 

p = -#(3p^ + 4p 7 ) . (2.11) 

It is convenient to work in terms of the dimensionless density parameters 



P P 
and the dimensionless "reduced" decay rate 



% = ^-, n 7 = ^, (2.12) 



° = vTh' (2 ' 13) 

which varies monotonically from to 1 in an expanding universe. 

The background equations (j2.9H2.ll]) can then be written as an autonomous system 

% = ^(n,-^-) , (2.14) 
n' = nJ-^—-n y ) , (2.15) 
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1(4 -0^(1-^ + ^(1 -<?), (2.16) 



where the prime denotes differentiation with respect to the number of e-foldings N = In a, 
and we have allowed for a time variation of the inflaton decay rate. The density parameters 
are subject to the constraint 

fy + fi 7 = l. (2.17) 
2.1 Gauge- invariant perturbations 

Linear scalar perturbations about a spatially-flat FRW background model are defined by 
the line element 

ds 2 = -(1 + 2>p)dt 2 + 2aB )i dtdx i + a 2 [(1 - 2ip)5 i:j + 2E tij ] dx i dx j , (2.18) 

where we have used the notation of Ref. [2T] for the gauge- dependent curvature perturba- 
tion, ip, the lapse function, ip, and scalar shear, x = a2 E — aB. 

The zero-th component of the perturbed energy transfer vectors, Eq. ()2.5|) . including 
terms up to first order, are written as [26] 

Q(<t>)o = -QiO- + V 2 ) - $Q<i> , 

Q( 7 )o = -Q 7 (l + <^)-5Q 7 , (2.19) 

where the gravitational redshift (time-dilation) factor (1 + ip) has been made manifest. 

Both the density perturbations, 5p ( f ) and £p 7 , and the curvature perturbation, if), are in 
general gauge-dependent. Specifically they depend upon the chosen time-slicing in an in- 
homogeneous universe. The curvature perturbation on fixed time hypersurfaces is a gauge- 
dependent quantity: after an arbitrary linear coordinate transformation, t — > t + 8t, it 
transfers as ip — > ip + H5t. For a scalar quantity, such as the energy density, the corre- 
sponding transformation is Sp — > 5p — pSt. However a gauge-invariant combination can 
be constructed which describes the density perturbation on uniform curvature slices or, 
equivalently the curvature of uniform density slices. 

The curvature perturbation on uniform total density hypersurfaces, (, is given by [27j 

( = -i>-H 5 4, (2.20) 
P 

while the curvature perturbation on uniform inflaton energy density and radiation energy 
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density hypersurfaces are respectively defined as 

P<t> 

C 7 = -ii-H 5 -^. (2.21) 

Pi 

The total curvature perturbation (|2.20|) is thus a weighted sum of the individual perturba- 
tions 

C = ^Cp + ^C 7 , (2.22) 
P P 

while the difference between the two curvature perturbations describes a relative gauge- 
invarariant entropy (or isocurvature) perturbation 

S H = m - Q = -3H ( 6 -P±- 5 -Pl) . (2.23) 

V P<f> Pi J 

From the definitions of the total curvature perturbation (|2.22j) and the entropy perturbation 
(|2.23|) . we get for instance that 

C* = C+^S H . (2.24) 

6 p 

On wavelengths larger than the horizon scale, the perturbed energy conservation equations 
for the inflaton energy density and the radiation energy density can be written, including 
energy transfer, as 



dp^ + 3H(5p (j} + 8P$) - (p^ + Py) 3ip = Q^ip + 

5p 1 + 3H(8 Pl + SP 1 )-(p 1 + P f )3ip = Q 1 ip + 8Q 1 . (2.25) 

The inflaton field and radiation have fixed equations of state (5P^ = and <5P 7 = 5p 1 /3) and 
hence there cannot be intrinsic non-adiabatic pressure perturbations. Using the perturbed 
(0 — ^-component of Einstein's equations for super-horizon wavelengths tp + Hip = — -f-^, 
we can re-write Eq. (J2.25J) in terms of the gauge-invariant curvature perturbations and 



7 Pi 
7 



where 



5Qintr,(f> = 8Qcf> —&P4> ) 

P<t> 

5Q intr , 7 = 5Q 7 -^5p 7 , (2.27) 
P 7 

are the gauge-invariant instrinsic non-adiabatic energy transfer perturbations, which are 
automatically vanishing if the local energy transfer functions Qi (i = <ft, 7) are functions of 
the local energy density p i; and 



^P \P<t> p) QHp 

SQrel,-y = ( 7~^" 7" J = ~ P^S^ (2.28) 



2p \ P 7 P / 6ifp 

are the gauge-invariant relative non-adiabatic energy transfer due to the presence of relative 
entropy perturbations 



2.2 Perturbing the inflaton decay rate 

Allowing for a perturbed decay rate (8T ^ 0) and for a possible time variation of T (T ^ 0), 
the perturbed energy transfer is simply given by 

5Q^ = -T5 P4> -5T P4> , (2.29) 
<5Q 7 = T5p^ + 5Y P(t> . (2.30) 

The corresponding intrinsic non-adiabatic energy transfer terms from the inflaton are given 

by 



<5Qintr, 7 = ^P<t> + r 5p0 - T — 5p4, - V— 5p 7 

p 7 p 7 



(2.3i; 
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The relative non-adiabatic energy transfer terms are given by 



SQmU = _r_M(Sj±_S_P\, (2 .32) 



2p \P<p p 
Tp^p ( <5p 7 dp 



<5Q re i, 7 = -f^ - -H . (2.33) 

2p V Pi P J 

Thus the evolution equations (|2.26|) for the curvature perturbation on uniform inflaton 
density hypersurfaces, and uniform radiation density hypersurfaces, £ 7 , are given by 

C, = —^S^ + H^STf, (2.34) 
6 P P4> P4> 



where 



3 p 7 V 2py p 7 



5Tf = ST-t 

P<t> 



ST? = 5T-r-^ = 5Tf + --^((-Q (2.36) 
7 Pj Hp 7 

are the gauge-invariant perturbations of the inflaton decay rate. As anticipated in the 
Introduction, the gauge-invariant decay rate perturbation receives contribution from two 
sources: the fluctuation of the decay rate 5T (generated during the inflationary stage) and 
the time variation of the decay rate T during reheating. 

Writing the total curvature perturbation on uniform total density hypersurfaces (J2.22)) 

as 

C = fC* + (1 - /) C 7 , (2.37) 

where 

_ 3f^ + gQ^ 

; " 4^ 7 + 3fi ' 1 4 * } 

we can obtain the following system governing the evolution of the total curvature pertur- 
bation on uniform total density hypersurfaces and the curvature perturbation on uniform 
inflaton density hypersurfaces 
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Q 4> (2g-3) r (i-g)n^ 



(C-C#) , 



(1 - g) (4 - fy) 4(1 -</)(!- fy) - ^ 




(2.39) 



In the standard inflationary scenario, p 7 is supposed to dominate the initial energy density 
after reheating and is assumed to be unperturbed at the end of inflation (beginning of the 
reheating stage) C 7jin = 0. Thus the curvature perturbation is initially (right after the end 
of inflation) an adiabatic density perturbation in the inflaton field. At this stage some 
comments are in order. 

• From the set of equations (|2.39J) . it is clear that during reheating, if 5T^ 1 = 0, the 
solution ( = is a fixed point attractor. Therefore, if at the end of inflation the 
total curvature perturbation ( is entirely provided by the curvature perturbation on 
uniform inflaton density hypersurfaces, d n = C<f>M, the total curvature perturbation 
( remains constant on super-horizon scales. At the end of the reheating stage and 
beginning of the radiation phase, C = C7 = (<j>M- 

• If the inflaton decay rate T depends only upon the inflaton field, T = T(0), then the 
gauge-invariant perturbation of the decay rate reads 



Since the long-wavelength solutions for the vacuum fluctuations in the inflaton field 



decay rate vanishes identically. This implies that the source term proportional to 
5T^ 1 for in Eqs. (I2.39J) vanishes identically and ( = is a fixed point attractor 

during the reheating stage. The total curvature perturbation ( remains constant on 
super-horizon scales during reheating. 

The same conclusion can be drawn if the decay rate is a function of the temperature 
associated to radiation. Suppose that the inflaton field is coupled to some fermion 
ip (radiation) through the Yukawa coupling Cy = hipipcf). This coupling allows the 




(2.40) 



obey the adiabatic condition 5<j)/<j) = 5<j)/<j), the gauge- invariant perturbation of the 
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inflaton field to decay into fermions with a decay rate T = ^M^y 1 — 4M|/M|, 

where and are the inflaton mass during the coherent oscillations and the 
fermion mass, respectively. During the reheating stage, at very early times, t <C 
r -1 , the energy density of the universe is dominated by the scalar field (f) and the 
radiation density is negligible. As the scalar field decays into fermions, the decay 
products rapidly thermalize forming a plasma with temperature T. The latter grows 
until it reaches a maximum value Tmax and then decreases as T oc a _3//8 up to the 
temperature Trh at the time t ~ T^ 1 which determines the end of reheating |2*9*1 I3U| . 
The thermalized fermions, produced during the first stages of reheating, acquire a 
plasma mass of the order of gT, where g is the typical (gauge) coupling governing 
the fermion interactions [SI]. This happens because forward scatterings of fermions 
do not change the distribution functions of particles, but modify their free dispersion 
relations, producing a plasma mass. The decay rate is therefore a function not only of 

the inflaton field 0, but also of the temperature T ~ plj A '. The gauge-invariant decay 
rate can be written as 



SI? = SI - m = fp 7 (l-'-?t) = (C - <*) , (2.41) 

which shows that the source term coming from ST^ 1 for ^ in Eqs. (12. 39(1 is propor- 
tioanl to (( — ^) and therefore ( = Q, remains a fixed point attractor during the 
reheating stage. The total curvature perturbation remains constant on super-horizon 
scales. 

If the inflaton decay rate depends upon the vacuum expectation value of another field 
X whose quantum fluctuations are excited during inflation and whose time variation 
during reheating is negligible, \ — 0, one recovers the inhomogeneous reheating sce- 
nario (TBI HE] with 5T^ 1 = 5T = (dT/dx) $X an d T = 0. During the reheating stage, 
when the inflaton field is still oscillating around the minimum of its potential and 
t r _1 , we can set ~ 1, g ~ 0. Eqs. (|2.39|) reduce to 

C * -(C-Q, 

C - - 577 , (2.42) 



11 



which is easily solved to give (going back to cosmic time) 



^^-\L M W;[/ fH{f)a{f) mf)' {2A3) 

where t in denotes the initial time of the reheating stage. Since on super-horizon scales 
we may consider 5T as a constant and a oc t 2//3 during the inflaton coherent oscillation 
phase, we obtain at the end of the reheating phase (which we set to be at t = 

Since deep in the radiation phase the gravitational potential is ^ 7 = — |C? we obtain 

A 5T 15T , A s 

ibry= ~ , 2.45 

^ 45 r 9 r ' K J 

which confirms the findings of Refs. fUJ QB] > m the case where the initial curvature 
perturbation ( in is tiny. 



2.3 Effects of a time variation of the inflaton decay rate 

Let us now analyze the new source of the variation of the total curvature perturbation on 
large scales proportional to T. We suppose that the decay rate is a function of time during 
the reheating stage, r ^ 0, and that its fluctuations during inflation are vanishing, 5T = 0. 
Under these circumstances, the gauge-invariant decay rate perturbation is given by 

5T Gi = _Y 6 -^. (2.46) 

Solving Eqs. (|2.39J) (going back to cosmic time), we find the total curvature perturbation 
at the end of the reheating stage is given by 

Cf = Cin- f dt' f(t')e J ^ f - 9)dt ' f dt"S(t")e- J i {f ' 9)dt " , (2.47) 

Jt in J tin 

where 



12 



./ 



H 



% (2g - 3) r (l-g)% 



(1 - g) (4 - H 2 4(1 -g)(l- fy) - ^ 
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g(4-fy) j 
2 (3 - 2g) 

P 1 - g <W 

3 + # /ty 



(2.48) 



Let us give an example. Suppose again that the inflaton field decays into some fermion ip 
(radiation) through the Yukawa coupling Ly = hipipcf). The fermion mass is a function 
of a scalar field S, 



The potential of the scalar field S reads V(S) = \M 2 S 2 + C 2 ^ (S - S f, with C> 1. 
If during inflation ^> Ms, the field S 1 is sitting at the position S = So and its quantum 
fluctuations are not excited since its effective mass ~ CH is much larger than the Hubble 
rate. At the end of inflation, the inflaton starts oscillating with mass M^. If XS > the 
inflaton cannot decay since the mass of the fermion is larger than M^/2 and the decay 
channel into fermions is kinematically forbidden. However, during the coherent oscillation 
phase, the Hubble parameter drops down as a -3 . If the Hubble rate becomes smaller than 
C~ 1 M$, the scalar field S may start rolling down towards the minimum at S = and settle 
there, thus decreasing the fermion mass. As soon as = XS becomes smaller than M^/2, 
say at some time t = t*, the inflaton decay channel into fermions becomes accessible and the 
decay rate rises from zero to a nonvanishing value Tq = ^M^. We can safely approximate 
the decay rate as T(t) — T 9 (t — i*) corresponding to 



From Eqs. (|2.39|) . we expect that the total curvature perturbation suffers a jump at t = £*. 
Inserting this expression in Eq. (J2.47)) and working in the limit Yq/H* < 1, we find that 
the total curvature perturbation jumps from the initial value Cin to the final value 



AL, = XS . 



(2.49) 



r 



T 5(t-Q . 



(2.50) 




(2.51) 
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In the opposite limit T /H* > 1, we find 



(2.52) 



One can envisage the alternative possibility that during the reheating stage the Hubble rate 
remains always larger than C^Ms- Under these circumstances, the field 5* does not roll 
towards 5 = 0, but remains stuck to the minimum of its potential at S = M ?^ H 2 So ■ The 

location of the minimum, however, changes adiabatically with time, S/S ~ —Q(Mg/C 2 H 2 ). 



AM 2 • X 2 S 2 M 2 

Correspondingly, the decay rate V = r W 1 — j^r changes with time as r/T ~ 48 ^^^ . 
Solving Eq. (|2.47|) gives 



Cf 





Sin — 







1 



96 \ 2 S 2 M 2 



15 C 2 M 2 Y 2 



Cn 



(2.53) 



We conclude that during the reheating stage the total curvature perturbation may be altered 
on super-horizon scales if the decay rate of the inflaton changes with time. Even though 
the amount of change is model dependent, the shift in the total curvature perturbation is 
always proportional to the value of the total curvature perturbation at the end of inflation 
and the beginning of the reheating stage Q n . 

2.3.1 Violation of the consistency relation 

During inflation both scalar and tensor perturbations are generated. As we have seen in the 
previous section, the total curvature perturbation can change on super-horizon scale during 
reheating if the inflaton field decay rate is a function of time. The difference between the 
values of the total curvature perturbation at the end of reheating and at the beginning of 
inflation can be parametrized by A = (Q — £i n ) / The power spectrum of scalar curvature 
perturbations at the end of the reheating stage and the beginning of the radiation phase is 
therefore given by 



V c (k) 



2tt 2 



Kf 



2tt 2 



;i + A) 2 Kini 2 



;i + a) 2 /H(k)Y 

e \Mp J 



(2.54) 
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where Mp is the Planck mass, H(k) indicates the value of the Hubble parameter when a 
given wavelength A = 2n/k crosses the horizon, i.e., when k = aH, and e = —H/H 2 is a 
slow-roll parameter accounting for the time variation of the Hubble rate during inflation. 

The primordial spectrum of gravitational waves is given by 



7T 2 V Mf t 

From expressions (|2.54|) and ()2.55|) . we can predict a consistency relation between the 
amplitude of the scalar perturbations, V^(k), the amplitude of the tensor perturbations, 
Vr(k), and the tensor spectral index, tit = dlnVT(k)/dlnk. Indeed, since Vx(k) oc H 2 (k), 
n T is given by n T = din H 2 {k)/ dink = — 2e and the consistency relation, for |A| <C 1, 
reads 

V T {k) 



V ( (k) 



16(l-2A)e= -8(l-2A)n T . (2.56) 



This consistency relation differs from the traditional one, Vt/V^ = — 8n^, obtained for 
one-single field models of inflation [32] when the inflaton decay rate during reheating does 
not change with time. Notice that, if A < 0, the tensor-to-scalar amplitude ratio can 
be larger than for one-single field models of inflation. Departures from the traditional 
consistency relation can be caused by other reasons: higher-order terms in the expansion in 
slow-roll parameters [22], quantum loop corrections [33] or the presence of multiple fields 
during inflation [33] . Given CMB B-mode polarization measurements, departures from the 
traditional consistency relation can be detected if the tensor-to-scalar amplitude ratio is 
larger than about 10~ 3 [33] and would rule out the simplest case of one-single field models 
of inflation with no variation of the decay rate. 



3 The nonperturbative inflaton decay 

As we already mentioned, if the classical inflaton field <fi very rapidly (explosively) decays 
into either its own quanta or into other bosons due to broad parametric resonance, the 
so-called preheating stage p!8j, Eqs. ()2.8|) should be modified. 

Suppose the inflaton is coupled to a light scalar \ with coupling C = \g 2 4> 2 "X 1 - During 
the coherent oscillations of the inflaton field, the x-q uan ta satisfy the so-called Mathieu 



15 



(n) , 

equation which leads to the existence of exponential instabilities Xk oc e^* * within the 
set of resonance bands of frequencies Ak^ labelled by an integer n and (for the first few 
bands) 



// 



(n) 



2 ' 



M q n 
2n {2 n ~\n- 1)!) 

For q^jl, preheating occurs mainly in the first resonance band with jjS 1 ' = and a 

width Ak^ ~ p 1 - 1 - 1 . For q ^> 1, however, many resonance bands are excited and preheating 

occurs in the broad resonant regime, with jj^ ~ 0.17, independent of q. At the beginning of 
the preheating stage, the so-called linear stage, the %-quanta grow exponentially in time till 
back-reaction effects set in. They are originated by the fact that the inflaton effective mass 
squared M 2 S = M| + g 2 (x 2 ) becomes dominated by the second term, thus decreasing the 
parameter q, and by the scattering of the produced particles off the zero mode (UHl E3 EE1 • 

During the linear stage, it is a good approximation to write the continuity equations 
(EZj) as 



P0 = -3H (p^ + P^) - pp 7 , 

p 7 = -3H (p 7 + P 7 ) + pp 7 , (3.2) 

where we have indicated by p 7 the energy density of light degrees of freedom, the x-particles 
generated during the first stage of preheating and by p the rate of production. Notice that 
p is a function of time as clear from Eqs. (|3.1|) . 

The background equations (j2.7J) can be re-written as 



% = fi 7 (fy - r) , (3.3) 
0; = rfi 7 fi , (3.4) 

r' = r^ + £(4-fy), (3.5) 
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where we have set 



r = § • (3.6) 

The local energy transfer functions are now given by = — pp 7 and Q 7 = pp 7 . Proceeding 
as in subsection 12. 1[ we can write the equations for the change of the curvature perturbation 
on uniform inflaton energy density and radiation energy density on super-horizon scales 



where 



P Pi ( 1 ^ Pi\ C l tjPi X..G1 



C* = ^(1-^)S H + H^f 1 (3.7) 
3 P<i> V 2 P / P4> 



C 7 =- ^ — (3.8) 
6 p p 7 p 7 



= 6/i- /i-^ , 



• ^7 _ r GI , A P 



^ = 5/i-/i-^ = 5/if + ^{C-Q (3.9) 
Pi H p 1 

are the gauge-invariant perturbations of the inflaton decay rate during the linear stage of 
preheating. 

The equations for the total curvature perturbation and the curvature perturbation on 
uniform inflaton energy density become 



, n , 3\ 3(^-r)(r-4) + |r(4-^)-^r[fi 7 (^-4)-l] 
C = <^( 2r -2) (r-4)(4-n,) 



From Eqs. IjH.lOj) . we conclude that if the resonance parameter p has either spatial or 
time variations, then the total curvature perturbation can be modified on large-scales. As 
a matter of fact, the parameter p is far from being constant during preheating. This 
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fact was not appreciated in previous studies of the effects of preheating on super-horizon 
perturbations (HH]- If the nonperturbative inflaton decay rate /J depends only upon the 
inflaton field, \i = /i(</>), then the gauge-invariant perturbation of the decay rate 5fi^ 1 
vanishes identically and ( = is a fixed point attractor during the preheating stage. The 
total curvature perturbation ( remains constant on super-horizon scales during preheating. 
This happens, for instance, if Eq. (jH.lj) holds and /x ~ q n ~ <p 2n . However, one can envisage 
other possibilities. For instance, as in the inhomogeneous reheating scenario, the coupling 
constant g 2 between the inflaton field <fi and the light bosons x can be a function of some 
other field which remains light during inflation. This leads to spatial fluctuations of the 
g-parameter. Tiny variations of q lead to large variations of the efficiency in extracting 
energy out of the inflaton zero mode [H7j and may have a dramatic effect on the final 
temperature fluctuations after the system thermalizes. This analysis will be presented in a 
separate publication [40J. 

4 Conclusions 

We have studied the evolution of large-scale curvature perturbations during the reheating 
stage after inflation in the case in which spatial and temporal variations of the inflaton 
decay rate are present. We have shown that the total curvature perturbation ( can change 
on large scales due to either spatial variations of the decay rate originated during inflation - 
the so-called inhomogeneous reheating scenario pHl ^] - or to a time variation of the decay 
rate during reheating. If only the latter source is present, the final curvature perturbation 
at the beginning of the radiation phase can be either smaller or larger than the curvature 
perturbation at the beginning of the reheating stage. This result leads to a violation of 
the consistency relation for one-single field models of inflation and to the observationally 
promising possibility that the tensor-to-scalar amplitude ratio is larger than in the standard 
scenario. Finally, we have presented the equations for the evolution of the total curvature 
perturbation and the curvature perturbation on uniform inflaton energy density hypersur- 
faces in the case in which the vacuum energy driving inflation is released into radiation 
through a (linear) stage of preheating. 
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